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ABSTRACT 

We present spectra of four Herbig AeBe stars obtained with the Infrared Spectrograph (IRS)^ on 
the Spitzer Space Telescope. All four of the sources show strong emission from polycyclic aromatic 
hydrocarbons (PAHs), with the 6.2 fim emission feature shifted to 6.3 /im and the strongest C— C 
skeletal- mode feature occuring at 7.9 /zm instead of at 7.7 /xm as is often seen. Remarkably, none of 
the four stars have silicate emission. The strength of the 7.9 fj,m feature varies with respect to the 
11.3 fiui feature among the sources, indicating that we have observed PAHs with a range of ionization 
fractions. The ionization fraction is higher for systems with hotter and brighter central stars. Two 
sources, HD 34282 and HD 169142, show emission features from aliphatic hydrocarbons at 6.85 and 
7.25 ^m. The spectrum of HD 141569 shows a previously undetected emission feature at 12.4 /xm 
which may be related to the 12.7 /im PAH feature. The spectrum of HD 135344, the coolest star 
in our sample, shows an unusual profile in the 7-9 /xm region, with the peak emission to the red of 
8.0 /xm and no 8.6 fim PAH feature. 

Subject headings: stars; chemically peculiar — infrared: stars 



I 1. INTRODUCTION 

Herbig AeBe stars (HAeBe) are intermediate mass 
(2-8 M(7)) analogs to the roughly so lar-mass T Tauri 

• stars plerbig.. „196nc iStrom et al.lll972D . The number of 
HAeBe stars in the Galaxy depends on the selection 
criteria, and as a result different studies have listed as 
few as ~50 (Herb ig's original list included 26 stars, but 
see Malfait ct al. 1998) or as many as several hundred 
ifrhe et al. 199 4'). One of the more conservative defini- 
tions assigns the HAeBe name to stars of spectral type A 
and B (and early F) that have broadened atomic emission 
lines in their optical spectra (hence the "e" in "AeBe") 
and have infrared emission in excess of a purely photo- 
sph eric spectral e nergy distribution. Using these crite- 

'. ria, iMalfait et alJ |l998) identified 45 HAeBe stars and 
pubhshed spectral energy distributions (SEDs) from UV 
to mm wavelengths, proposing an evolutionary scenario 
base d on mid-i n frared colors and SED structure. 

Meeus et al.l l|2001|) observed the 14 spatially unre- 
solved sources (at all wavelengths) from the list by 
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IMalfait et all l|1998t) with the Short- Wavelength Spec- 
trometer (SWS) aboard the Infrared Space Observatory 
(ISO). They report the detection of polycyclic aromatic 
hydrocarbons (PAHs) in seven of the 14 sources, propos- 
ing that the emitting PAHs are located in flared disks 
where the shadow of a t hicker inner disk cannot atten- 
uate the UV radiation. lAcke fc van den Anckeil l|2004f) 
recently presente d spectra f rom the SWS for the entire 
list of sources of IMalfait et al. (1998), and they report 
the detection of PAH emission in many of them, men- 
tioning the possibility that differences in the ionization 
fraction of the PAHs could help explain the structure 
in the 2-20 /im spectra. They also noted that emission 
from crystalline silicates may confuse the issue. It is also 
clear from the SWS spectr um of 51 Oph (HD 156843; 
Ivan den Ancker et alJl20?)ll) that emission from gas near 
the stars may play a significant role in shaping the mid- 
infrared spectra of these stars and complicating detailed 
st udies of PAH emissi on in that spectral range. 

iHabart et alJ l)2004j ) have analyzed ground-based and 
ISO observations of some 30 HAeBe stars and applied a 
model in which the material in the disk is in vertical hy- 
drostatic equilibrium, is heated by the central star, and 
contains large dust grains in thermal equilibrium with 
the radiation field as well as small grains and PAHs that 
are transiently heated. They suggest that the PAH emis- 
sion is primarily from the outer disks (R~100 AU) and 
that strong PAH emission is predominantly from sources 
that have fiared disks (as earlier proposed bv lMeeus et alJ 
•Ml*). 

Li fc Luninj l)2003j) modeled ground-based spectra of 
HD 141569 in the 8-13 and 17-25 /xm regions using a 
porous, cometary-type dust model consisting of coagu- 
lated, but otherwise unaltered interstellar grains along 
with fully ionized PAHs. They predict the shape of the 
5-25 /im spectrum, noting that there should be little 
emission from crystalline silicates (no more than 10% 
by mass fraction). Since PAHs are a significant compo- 
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nent of the interstell ar medium (ISM ), important to the 
energy balance fe.g. IAllamandola~et al. 198^, and are 
apparently necessary to accurately model the emission 
in HAeBe stars, it is time to closely investigate not just 
their presence, but the specific spectral properties of the 
emission features from PAHs around HAeBe stars. 

High-sensitivity spectra of PAHs allowing precise band 
strength measurements can provide information about 
the structure of disks and their physical and chemical 
conditions. The presence or absence of PAH features 
may provide clues to the physical properties of the disks 
in HAeBe systern s, such as their degree of flaring (cf. 
iMeeus et al.ll200l|) . and may also serve as an evolution- 
ary clock that tracks grain processing, since there are 
multiple destruction paths for PAHs in disks (intense ir- 
radiation, compaction and incorporation in larger grains) 
but no clear formation path. The PAHs in disks around 
HAeBe stars most likely originate in the ISM from which 
the stars formed. As a significant interstellar component, 
PAHs may have condensed onto the ice mantles of dust 
grains in dense clouds that were then incorporated into 
protostellar nebulae. The PAHs emerge as free-flying 
molecules when the ice mantles of grains sub limate due to 
stella r irradiation or grain-grain collisions ijLi fc Lunind 
12001 . The origin of interstellar PAHs, however, is not 
well constrained. Suggested sources for interstellar PAHs 
includ e: (1) format ion around and ejection from carbon 
stars ()Latteilll99H) . (2) shattering of carbonaceous in- 
terstellar dust, or of photoprocessing of man tles on in- 
terstellar dust grains (arce nberg et a l.'2000^. by gram- 
grain collisions in interstellar shocks (Jones ct al. 1996), 
or (3) in-situ formation through ion-molecule reactions 
(|Herbsdl199T]) . It is essential to quantify the PAH fea- 
tures as a diagnostic in order to distinguish between these 
possibilities. 

We have begun a spectroscopic survey in the mid- 
infrared (5.2-3 6 um) of 16 HAeBe stars from the sam- 
ple defined by iMalfait et all lll998D using the Infrared 
Spectrograph (I RS: iHouck et all I2004D on the Spitzer 
Space Telescope l)Werner ^t'lLLfl2004t) . One of the goals 
of this study is to closely examine the PAH emission 
in the 5.2-14 /im spectral region. We present the first 
results of our study with spectra and analysis for four 
stars that show PAH features, no silicate features, and 
few other solid state features in their spectra. This has 
allowed a careful examination of the relative positions, 
strengths, and shapes of the features at 6.2, 7.7-7.9, 
11.3, and 12.7 fj,m without concern for the effects of sili- 
cate dust grains, either amorphous or crystalline. While 
ISO/SWS and PHOT-S spectra of the four sources in our 
sample show PAH emission fe atures (Mccus ct al. 2001; 
lAcke fc van den Anckeill2004(l . the Spitzer IRS spectra 
have higher sensitivity, allowing the detailed analysis 
that we present here. 

All four of our sources sh ow evidence of circumstellar 
disks in their SEDs (Malfait 'eral]ll998|) . HD 34282 evi- 
dently has a large Keplerian disk inf erred from millime - 
ter observations l|Merin et al. 2004: IPietu et alJl200l . 
and the disk around HD 141569 has been imaged di- 
rectlv in the optical and near infrar ed llBoccate li et al.l 
1 20031 iMouillet et a l. 2001; Wei nberger et all I1999F 
lAugereau et al.l 119991 Weinber ger et all '2005). HP 
135344 ha s yet to be s p atially resolved (Mcc us et al.l 
120011) . hut iHabart et alJ l|2005ri have spatially resolved 



the PAH emission at 3.3 ^m in HD 169142. 

2. OBSERVATIONS AND ANALYSIS 

The four HAeBe targets were observed with the IRS 
on Spitzer in 2004 from January to October, as described 
in Table 1. This paper focuses on the observations with 
the Short-Low module (SL), which covers the 5.2-14 fim 
region with a spectral resolving power of R~90. The SL 
slit is 3'.'6 (2 pixels) across and 57" long. 

We started with the flatfielded images generated by 
the standard IRS data reduction pipeline (Sll.O) at the 
Spitzer Science Center (SSC). We removed the back- 
ground emission by differencing images with the source 
in the separate nod positions of each of the SL aper- 
tures. When extracting spectra from the images, we used 
an extraction window which scales with wavelength; this 
method is similar to that used by the SSC pipehne. We 
obtained a spectrophotometric calibration by using iden- 
tically extracted spectra of the standard star HR 6348 
(KO HI) and spectral templates of t his star p repared us- 
ing the methods descri bed bv iCo hen et al,| (12003) and 
modified as described bv lSloan et al. (2005b). The spec- 
tra in both nods and in both SL apertures were calibrated 
separately and then combined into one spectrum. This 
method generally results in very high precision within a 
single spectrum, indicated by the error bars in our plots, 
and an absolute photometric accuracy of ^--^10% or better. 

We have modified the Sll. O version of t he wav elength 
calibration, as described by iSloan e~ D ll2005aD . This 
modification shifts the wavelength calibration in SL or- 
der 1 by 0.04 /im to the blue. We estimate that our 
wavelengths are accurate to better than 0.02 /im. 

Figure 1 presents the spectra of our four HAeBe stars 
(with the spectrum of HD 135344 divided by 10 to fit on 
the same plot with the other three). The figure includes 
error bars which are generally smaller than the vertical 
widths of the lines used to plot the spectra. These errors 
are the formal uncertainty in the mean of the spectra 
from the two nod positions in the SL slit (half the ab- 
solute value of the difference). With the notable excep- 
tion of HD 141569, which appears to have a single cool 
component peaking at 40 /^m, the SEDs in our sample 
consist of a combination of a warm component peaking 
in the near infrared (^2 fim) and a cooler comp onent 
peaking at longer wavelengths (IMalf ait et al.ll9 9ffl. This 
double-peaked structure may indicate that the disks have 
formed significant gaps between a warm inner disk and 
a cool outer disk. The wavelength of the dip in intensity 
between the warm and cool components shifts slightly 
from star to star. Thus three of the spectra are flat be- 
cause the wavelength region observed by SL lies more in 
the dip between these two components. The continuum 
of HD 135344 appears blue because its dip appears to- 
wards longer wavelengths than the other sources and is 
therefore dominated by the warm component (T ^ 1000 
K) at the short wavelength end of our spectrum. 

Figure 2 illustrates how we measure the strengths of 
the PAH features, using HD 141569 as an example. We 
have used line segments as a rough estimate of the con- 
tinuum under each feature. We anchor the line segment 
to the spectrum on either side of each feature and in- 
tegrate the residual spectrum between the two anchor 
points to determine the total flux in the feature (in W 
m^^). Table 2 lists the wavelength intervals used to fit 
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TABLE 1 
Observing information 



Target 


Spectral Type 


U{r) (x 10^)=" 


IRS Campaign 


Observation Date (UT) 


AOR Key 


HD 34282 


AOe 


2.3 


13 


2004 Oct. 10 


3577856 


HD 141569 


AO Vei^ 


2.3 


4 


2004 Mar. 3 


3560960 


HD 169142 


A5 Ve'' 


0.71 


5 


2004 Mar. 26 


3587584 


HD 135344 


F4 Ve*^ 


0.16 


11 


2004 Aug. 8 


3580672 



'^U{r) is the radiation field due to the star in units of the interstellar radiation field in the Solar neighborhood, calculated at r=100 AU 
''Revised spectral type from lDunkin et al.l 119971) . 
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Fig. 1. — Spectra of the sample of Herbig AeBe stars showing 
PAH emission features and no features from silicate dust grains in 
their spectra. The spectrum of HD 135344 has been divided by 10. 
These spectra do include error bars, although they are generally 
smaller than the width of the lines used to plot the data. 
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Fig. 2. — The spectrum of HD 141569 showing the linear fitting 
method used to measure the equivalent flux in the PAH features. 
Table 2 gives the wavelengths used to anchor the line segments on 
either side of the PAH features. 



each feature. This method assumes that the PAH emis- 
sion is optically thin, and it avoids dependencies on more 
sophisticated techniques to model the continuum. HD 



TABLE 2 
Fitting wavelengths 



Feature {fJ.m) 




Xred (/im) 


6.2 


5.60-5.90 


6.71-6.78 


7.9 


6.98-7.07 


9.20-9.40" 


11.3 


10.83-10.90 


11.70-11.86 


12.7 


12.10-12.30 


12.97-13.10 



=^8.85-9.05 fim for HD 135344. 



135344 is the exception; before extracting the features 
using line segments, we first fit a polynomial estimate of 
the continuum to the spectrum between the features and 
subtract it. 

We have avoided the problem of separating the PAH 
features from other dust features like amorphous sili- 
cate emission at 10 fim or crystalline silicate emission 
at 11.2 /im by choosing our sample carefully. The pri- 
mary uncertainty is how much of the emission below the 
line segments belongs to the underlying PAH emission 
plateaux and how much belongs to the feature we are 
measuring. If we are truncating flux from a PAH fea- 
ture, we are doing so consistently to all of the stars in 
the sample, and our relative measurements should re- 
main valid. 

The resulting profiles for the PAH features at 6.2, 7.7- 
7.9, 11.3, and 12.7 /xm appear in Figures 3 and 4. Table 
3 presents the total fluxes and central wavelengths for 
each PAH feature in each of our sources. The nominal 
7.7-7.9 /im feature is shifted to 7.9 fim (or further) in our 
sample, so for the remainder of the paper, we will refer 
to this feature as the 7.9 /im feature. The reader should 
note that the flux labelled Fy.g in Table 3 and plotted in 
following figures includes the emission of the 8.6 /im fea- 
ture as well. It could be more properly described as the 
total PAH flux from 7 to 9 /zm. We estimated the central 
wavelength of each feature by finding the wavelength at 
which half the emission lies to the blue and the red and 
by fitting Gaussians. The two methods produce similar 
results, but we adopt the Gaussian fits because they are 
less susceptible to noise and blended features. 

Since ultraviolet and visible photons fr om the stel- 
lar p hotosphere excite the PAH emission ijLi fc Draind 
120021) ■ we have estimated integrated fluxes for each of 
our program stars, relative to the interstella r radiation 
field (ISRF) in the Solar neighborhood (,,Mathis. et alJ 
11983(1 . These estimates are intended only as a guide in 
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TABLE 3 
Feature positions and strengths 





f 


3.2 /im feature 




7.9 fim feature 


11.3 


feature 


12.7 


fim feature 




A 


Flux 


A 


Flux 


A 


Flux 


A 


Flux 


Target 


(fim) 


(10-15 W m-2) 


(/im) 


(10-15 W m-2) 


(/im) 


(10-15 W m-2) 


ifim) 


(10-15 w m-2) 


HD 34282 


6.26 


12.63 ±0.17 


7.92 


47.09 ±0.58 


11.27 


3.11 ±0.14 


12.64 


0.54 ±0.06 


HD 141569 


6.26 


6.84 ±0.15 


7.98 


16.65 ±0.35 


11.28 


0.62 ±0.10 


12.52 


0.18 ±0.02 


HD 169142 


6.26 


32.71 ±0.54 


8.01 


88.06 ±0.93 


11.30 


11.86 ±0.58 


12.71 


1.95 ±0.17 


HD 135344 


6.28 


3.47 ±0.51 


8.02 


12.49 ±1.30 


11.36 


1.38 ±0.21 


12.75 


0.18 ±0.06 




A (um) 

Fig. 3. — The 6.2 fim PAH emission feature and the 7-9 
PAH emission complex (with error bars), all normalized to the 
peak emission in the vicinity of 8 ^m. Line segments have been fit 
under the two emission regions and subtracted. In all four spectra, 
the "6.2 fira" feature has shifted to 6.3 fira, and the peak emission 
around 8 /im occurs at (or beyond) 7.9 fim and not 7.7 ^m as 
often seen in other objects. The spectrum of HD 135344 shows 
an emission plateau from 8.0 to 8.2 fim and is missing the 8.6 fim 
feature. 
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Fig. 4. — The PAH emission features at 11.3 and 12.7 fim (with 
error bars), after continuum subtraction. Separate line segments 
have been fit under the features and removed. All spectra are 
normalized at 11.3 fim. The most ionized spectrum, HD 141569, 
shows an emission feature at 12.4 fim. 



our discussions of PAH excitation. iLi fc Lunind 1)200.'^ 
calculated the integrated flux (A=0.09-1.00 ^m) for HD 
141569 as a function of distance, r, from the central star 
using spectral type B9.5 V and a Kurucz photosphcric 
model. Table 1 presents the integrated fluxes of our 
stars for r=100 AU calculated relative to their value for 



HD 14 1569. Following the nomenclature of ILi fc Lunind 
lj2003ll . U{r) = 1 corresponds to the ISRF in the solar 
neighborhood. We note that these are not in units of 
Habings or Gq, the ratio of the UV energy density to the 
estimate of Habing iHabing 1 119681 ) . since we include the 
visible spectrum in addition to UV. Because the PAHs 
are emitting either from the surface of the disk or from 
the surrounding envelope, we assume no absorption be- 
tween the photosphere and the PAH molecules. 

3. DISCUSSION 
3.1. General spectral characteristics 

The 7.7 fim PAH feature consists of two main com- 
ponents, one at 7.65 fim a nd the other at 7. 85 fim 
ijCohen et al1 IT989: Bregma^ [T989t iPeeters erali r2002). 
The 7.65 /im feature tends to dominate the PAH spectra 
of reflection nebulae and H II regions where the PAHs 
seem to have been heavily processed. The 7.85 fim fea- 
ture dominates in many planetary nebulae and objects 
evolving away from the asymptotic giant branch (AGB) 
where the PAHs seem to be relatively fresh and unpro- 
cessed, peters et al. (2002) studied 57 spectra of PAH 
sources obtained with the SWS, and they found that in 
spectra where the 7.85 /im component dominated, the 
6.2 fim PAH feature was shifted to 6.3 fim. They used 
these two characteristics as the basis for what they called 
Class B PAH spectra. 

Of the 57 spectra studied bv IPeeters et alJ 1)20021) . 42 
belong to Class A, which shows the PAH emission fea- 
tures at 6.2 and 7.7 /im. Class B includes 12 spectra, and 
Class C, where the peak emission in the 7-9 /zm range is 
shifted to the red of 8.0 /im, accounts for 2 sources. One 
source in their sample had characteristics of both Class 
A and B. 

Figure 5 compares the spectra of HD 141569 and HD 
135344 in our sample to some of the sources considered 
by Pcetcrs ct al. (2002). The top panel illustrates PAH 
spectr a of Classes A, B, and C as deflned bv lPeeters et alJ 
1)2002(1 ■ using as prototypes IRAS 03260-f3111 (a young 
stellar object in the NGC 1333 molecular cloud core), 
HD 44719, (the Red Rectangle), and AFGL 2688 (the 
Cygnus Egg) , respectively. Th e comparison spect ra were 
obtained from the SWS Atlas l)Sloan et al .' 2003). AU of 
the IRS spectra and the SWS spectra of IRAS 03260 and 
HD 44179 have had the continuum removed using line 
segments as illustrated in Figure 2. The SWS spectra of 
AFGL 2688 and HD 100546 have steep red continua; to 
isolate the PAH spectrum we flrst fit a polynomial to the 
continuum and subtracted it. 

Figure 3 and the middle panel of Figure 5 show that 



PAHs in Herbig AeBe Stars 



5 




6.0 6.5 7.0 7.5 8.0 8.5 9.0 

Fig. 5. — A comparison of spectra from the SWS on ISO with two 
spectra from our samp le. Top: P rototypical SWS spectra of the 
three classes defined bv lPeeters e t al. 12002). IRAS 03260+3111 is 
a young stellar object, while HD 44179 (the Red Rectangle) and 
AFGL 2688 (the Cygnus Egg) are post-AGB objects. Middle: HD 
141569 compared to two Class B objects. HD 141569 resembles the 
HAeBe star more than HD 44179. Bottom: HD 135344 compared 
to HD 100546 and the Class C spectrum from AFGL 2688, showing 
that the spectrum resembles a blend between the two. 



our HAeBe stars have spectral features typical of Class 
B, with the 6.2 fj,m feature shifted to at least 6.26 /im, 
and the peak emission in the 7.6-7.9 /zm region falling 
near 7.9 fim. Both the 6.2 and 7.9 /Ltm features arise 
from C— C stretching modes. The 8.6 /zm feature, which 
is included in our measurements of the strength of the 
7.9 fim feature, probably arises from a C— H in-plane 
bending mode. 

The sample examined bv lPeeters et al.l l|2002|) includes 
two isolated HAeBe stars, HD 100546 and HD 179218, 
both of which they included with the Class B PAH spec- 
tra. The middle panel of Figure 5 shows that in the spec- 
trum of HD 100546, the 7.9 ^m PAH feature is actually 
shifted to even longer wavelengths, nearly to 8.0 fim. In 
the other HAeBe star, this feature appears at 7.8 /im, 
but in all four of the IRS spectra, this feature lies closer 
to 8.0 /im. Thus, i n five of the s i x isola ted HAeBe stars 
studied here or bv lPeeters et alJ l|2002|) . the 7.8-7.9 fiui 
feature indicative of Class B is shifted nearly to 8.0 /zm. 
The HAeBe stars may represent a subgroup within Class 
B. 

Other features visible in the spectra in Figure 1 in the 
5-8 /im range require comment as well. The 5.2 /im PAH 
emission feature appears in three of the spectra at the 



short-wavelength cut-off, but it is absent in HD 34282. 
Two of the spectra, HD 34282 and HD 169142, show 
spectral structure in the vicinity of 7 /im which is due to 
em ission features at 6.85 and 7.25 /im. 

Chi ar et al.l ((2000) made the first detection of absorp- 
tion features at 6.85 and 7.25 /im, which they attributed 
to C— H stretching modes in aliphatic hydrocarbons such 
as methyl or methylene groups. This discovery was made 
in the SWS spectrum of Sgr A*. Bouwman et al. (200]]) 
detected the 6.85 /xm ban d in emission in the HAeBe star 
HD 163296. iSooon et al. (2004) recently reported the 
6.85 and 7.25 /xm bands in absorption in the spectrum 
of the ultraluminous infrared galaxy IRAS F00183-7111, 
and they showed that lab oratory data of hyd rogenated 
amorphous carbon (HACs; iFurton et al.l ll999^ provide a 
good match. We note that the 3.4 /im interstellar ab- 
sorpti on band, identified as an aliphatic C— H stretching 
mode ("Sa ndford et alJll991ll. can also appear in emission 
(Ccballe fcTan der Ween lll^ iGebaile et alJll99l . It 
is likely, therefore, that we have observed these aliphatic 
HAC features in emission in HD 34282 and HD 169142. 



3.2. HD 135344 

The spectrum of HD 135344 contains a much stronger 
contribution from a warm (T ~ 1000 K) underlying the 
PAHs, making the shape and position of the weaker PAH 
features much more susceptible to any systematic prob- 
lems with the extraction and calibration of the spectrum, 
which explains the larger error bars in its spectrum in 
Figure 3. Nonetheless, some differences from the other 
spectra in the sample are still apparent. 

In the 7-9 /im region, HD 135344 shows a plateau ex- 
tending from 8.0 to 8.2 /im, most likely a blend of peaks 
at these two wavelengths. The 8.0 /im peak is similar 
to that seen in four of the other HAeBe stars, but the 
8.2 /im peak more closely resembles Class C PAH spectra 
(^eetcrs et al. 2002). Both of the objects in this class, 
AFGL 2688 and IRAS 13416-6243, are post-AGB ob- 
jects, and neither show an 8.6 /im feature. HD 135344 
has more emission shortward of 8.0 /im than these two 
objects, but the missing 8.6 /im feature is striking. 

The spectrum of HD 135344 differs from the other 
HAeBe stars in other ways. The 11.3 /im feature is 
shifted by at least 0.06 /im to the red; this is too large 
to arise from the random pointing errors which occasion- 
ally place a star off the central axis of the IRS slit. HD 
135344 also has a relatively weak contribution 6.0 /im 
(compared to 6.3 /im). 

Beintima et al. (1996) attributed the 6.0 /tm satellite 
of the 6.2 /im feature to a weak C— C mode, but we note 
that in the objects studied bv iPeeters et al.l l|2002(l the 
strength of the 6.0 /im feature shows no correlation with 
that of the 6.2 /im C— C mode. A correlation between 
the stre ngth of the two ban ds also seems absent in our 
sample. IPeeters et al.l (|2002j) speculated that the 6.0 /im 
feature arises from C— O stretching modes in oxygenated 
PAH species, but these modes generally peak at some- 
what shorter wavelengths l|Schutte et alJll990|) . The car- 
ri er of the 6.0 /im feature remains uncertain. 

lOudmaiier et al.l |1992^ have revised the spectral clas- 
sification of HD 13534 4 from AO V (,Houk 198a) to F4 
Ve, and iDunkin et al.l l|1997() have confirmed this new 
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spectral type.^^ It is interesting that this source, which 
is significantly cooler than the other three sources in the 
sample, is the one which shows the most unusual PAH 
emission in the 8 ^m region. 

3.3. Dependence of PAH emission on ionization fraction 

The relative strengths of PAH features can be used to 
identify the degree of ionization of the PAHs. They may 
therefore provide clues to the incident radiation field, 
electron densisty, and gas temperature. The most ob- 
vious difference between neutral and ionized PAHs is the 
relative strength of the C— C modes in the 7.7-7.9 /im re- 
gion to the out-of-plane C— H bending modes at 11.3 /im 
and longer wavelengths (e.g. lAilamandola ct al, ,199 9) . 
Ionized PAHs emit more strongly in the features at 6- 
9 /im than in the 10-13 /Ltm region, while in neutral 
PAHs, the 1 -13 ju m region dominates. The models of 
iDraine fc Lil l)200Hl show that the ratio Fy.g/i^n.a in- 
creases with the size of the PAH molecule since larger 
PAHs generally h ave a higher C:H ratio (e.g. see 
iLi fc Drainell200'l|l . However, a shift from neutrals to 
ions increases the ratio much more significantly, by a 
factor of ~10. 

We can thus use the ratio of the strengths of the 7.9 
and 11.3 /^m features as an estimate for the relative de- 
gree of ionization in each spectrum. This ratio is the 
most obvious difference among the spectra in our sam- 
ple, varying from 7.4 in HD 169142 to 25.2 in HD 141569, 
making these the least and most ionized PAH spectra, re- 
spectively. The ratios for the intermediate stars are 9.1 
for HD 135344 and 15.2 for HD 34282^ 

Recent optical spectroscopy by iDunkin et al.l l)1997fl 
confirms the previous revision of the spectral classifica- 
tion of HD 135344, as discussed above, and also revises 
the classification of HD 169142 from B9 V ( Houk 1982) 
to A5 Ve. These new classifications correlate with the ob- 
served ionization fractions as measured by the -F7.9/F11.3 
ratio, with the two AO stars having the more ionized 
spectra, and the later-type stars having the less ionized 
spectra. 

In HD 1 41569, according to detailed models 
ijLi fc Luninel l2003). the charging of PAHs is dominated 
by collisions with electrons and not by emitting photo- 
electrons. The electron recombination rate is three or- 
ders of magnitude larger than the photoelectron emis- 
sion rate. For this case, then, the PAHs are negatively 
charged. Our analysis is not sensitive to the polarity of 
the charge of the PAHs, only the ionization fraction. 

3.4. The 7.9 fim feature 

There is no clear explanation for the varying strengths 
of the components of the 7-9 jum PA H complex among 
observed sources. iBreeman fc Teml (j2005l|) investigate 
two possible causes in their study of the spatial distri- 
bution of the components in extended reflection nebu- 
lae. The peak of the complex might shift from 7.6 to 

11 IZuckcrman ct alj <T995h and ICoulson fc Waltheil 
find ev en later spec tral types, F6 V and F8 V, respectively. 
ICoulson fc Waitheii jim^ point out that HD 135344 is a visual 
binary with the two components, SAO 206462 and 206463, sepa- 
rated by 20". SAO 206463 is the brighter source and is probably 
the origin of the AO V classification, but the coordinates of HD 
135344 match those of the fainter SAO 206462, which is the target 
under consideration here. 



7.9 /im as the fraction of PAH cations increases or as 
the fraction of fresh and unprocessed PAHs increases. 
Our sample shows a range of ionization fractions, but in 
three of the four spectra, the shape of the 7-9 fj,m com- 
plex looks similar. The ratio of the 7.9 /^m and 11.3 /im 
features varies by a factor of 3.7 in these three sources, 
and yet the relative contributions of the 7.6 and 7.9 iim 
components seems fixed, which would suggest that the 
components are not sensitive to ionization, at least over 
the ionization range we have examined. 

The size of the PAH molecules may also influence the 
relative strengths of the features at 7.6 and 7.9 /im, but 
there is some disagreement over how the position of the 
strongest featur e in this wav e length range will depend 
on PAH size. iBoissel et al.l l|2001fl suggest that very 
small grains produce t he longer- wavelength feature, but 
iHudgins et al.l l)1999b|) report that their laboratory mea- 
surements show that the 7.6 /im feature shifts to longer 
wavelength for larger PAHs. We can only say that our 
data are least consistent with the dependence of the 7- 
9 /tm PAH complex on ionization fraction. 

3.5. Possible correlations with PAH size and structure 

The models of IDraine fc Lil l|2001f) show a dependence 
of the ratio F6.2/-F7.9 on the size of the PAH molecule, 
but it is weak, with the flux ratio doubling only when the 
number of C atoms has been reduced by a factor of ~10. 
The top panel of Figure 6 shows that the ratio ^6.2/^7.9 
varies by a factor of two in our sample, but not in any 
systematic way with the ionization fraction as indicated 
by F7.9/F11.3. 

To better discri minate the size of the PAH molecules, 
iHonv et all l)200H) use the ratio Fi 9 7 /Fi 1 The 11.3 /tm 
feature arises from a C— H stretching mode when an aro- 
matic ring contains only one H atom, such as one might 
find along the edge of a large PAH molecule. The 12.7 /tm 
feature arises from the C— H trio stretching mode, i.e. 
from aromatic rings with three H atoms, which will be 
found on the corners and ends of PAH molecules. As 
one increases the size of a PAH molecule, the strength of 
the emission at 11.3 /im will increase with respect to the 
emission at 12.7 /tm. 

The bottom panel of Figure 6 shows that the strength 
of the 12.7 /im feature generally increases with the ion- 
ization fraction. Comparing this panel with the top 
panel reveals that the use of the ratios ^6.2/^7.9 and 
^i2.7/-Fii.3 as size indicators gives contradictory results. 

The -^12. 7/^11. 3 ratio may be influenced not just by the 
size of the PAHs, but by their geometry. If a compact 
PAH is rearranged into a structure with the same num- 
ber of carbon atoms but with more uneven edges, the 
strength of the duo and trio modes will increase com- 
pared to the solo mode. Thus, it is possible that the the 
enhancement of the 12.7 /tm feature could also arise from 
more irregularly shaped and less compact PAHs. 

Measurements at the Astrochemistry Laboratory at 
NASA Ames Research Center have revealed no depen- 
dency of the ^6.2/^7.9 ratio on the size of the PAHs (Hud- 
gins 2004, private communication), but the Astrochem- 
istry Laboratory has not explored a full order of mag- 
nitude in P AH size. While lar ger PAHs, such as those 
modeled bv D raine fc Lil lj2001(l . do show a dependence, 
typi cal PAH sizes are o n the order of ^100 carbon atoms 
(e.g. lLi fc Drainell2001(l . and the range in PAH sizes may 
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Fig. 6. — Flux ratios from the different PAH features plotted vs. 
the ratio of the 7.9 and 11.3 /im PAH features. The ratio -Fr.g/i^ii.s 
traces the ionization fraction of the PAHs. From left to right (least 
to most ionized), the symbols are for HD 169142 (A5), HD 135344 
(F4), HD 34282 (AO), and HD 141569 (AO). In the top panel, the 
relative strengths of the 6.2 and 7.9 fim features vary by a factor of 
~2, but they do not vary systematically with ionization fraction. 
The bottom panel shows that the ratio F12.7 / Fu s increases with 
ionization fraction, suggesting that the size of the PAHs decreases 
with ionization fraction. 

be insufficient to dominate other causes for variation in 
tfie strength of the 6.2 and 7.9 fim bands. The relative 
strengths of the out-of-plane C— H bending modes, on 
the other hand, are much more sensitive to variations in 
PAH size, especially for the sizes of the molecules likely 
to be encountered in typical astrophysical environments. 

The spectra of our HAeBe stars show other emission 
features in the 11-13 /im region in addition to the 11.3 
and 12.7 /xm features. In Figures 1 and 4, the duet 
out-of-plane C— H bending mode at 11.9 /im is visible 
in all of the spectra except HD 135344. These figures 
also show that what we are measuring as the 12.7 /Ltm 
feature shows some structure. In the spectrum of HD 
141569, the secondary emission peak at 12.4 /im actually 
dominates the 12.7 fim feature. The Humphreys-a line 
should be at 12.365 /im, but there is no corresponding 
Pfund-a hue at 7.46 /im. The H2 0-0 S(2) is expected 
at 12.28 /im, but this is too far to the blue and it is 
not accompanied by the 0-0 S(3) line at 9.67 /im. The 
most likely carrier of the 12.4 /im feature is another trio 
out-of-plane bending mode, as this feature can vary from 
12.1 to 13.5 /im, depending o n the PAH molecule (e.g. 
iHudgins fc Allamandolall 1999(1 . 



4. CONCLUSION 

We have presented mid-infrared spectroscopy from 
Spitzer of four Herbig AeBe stars whose spectra show 
strong PAH emission but no silicate dust features. The 
absence of silicate emission has enabled a detailed anal- 
ysis independent of specific disk models of the relative 
strengths and positions of PAH emission features in the 
5-14 /im region. 

Our sample shows spectral characteristics generally 
consist ent with the Class B PAH emission sources identi- 
fied bv lPeeters et~al] l|200l . viz. a 6.2 /im feature shifted 
to 6.3 /im and 7-9 /im emission dominated by the com- 
ponent at 7.9 /im. However, in all of our sources, the 
7.9 /im feature is shifted closer to 8.0 /im, as in the spec- 
trum of HD 100546, a HAeBe star observed by the SWS 
on ISO. 

The spectra show a range in ionization fraction as re- 
vealed by the i^r.g/Fn.s ratio, with HD 141569 showing 
the highest degree of ionization and HD 169142 the low- 
est. The ionization fraction is lower for later spectral 
classes. The strength of the 12.7 /im feature increases 
with respect to the 11.3 /im feature as the ionization 
fraction increases. The strength of the 6.2 /im feature 
relative to the 7.9 /im feature varies by a factor of nearly 
two in our sample, but in no systematic way. 

The spectrum of HD 135344, the coolest star in our 
sample, resenibles the Class C spectra discussed by 
iPeeters et all (|2002^ in some ways, with more emission 
to the red of 7.9 /im than in the other sources and a 
missing 8.6 /im feature. 

Thus, our sample of PAH spectra from HAeBe stars 
show some characteristics which reinforce results from 
the literature, but also raise new questions. Why do 
HAeBe stars consistently produce the Class B PAH spec- 
tra seen in less than 25% of the SWS sample? What 
causes the shift of the 7.9 /im feature to ~8.0 /im? What 
causes the range of 6.2 /im feature strengths compared 
to the 7.9 /im feature? What is the carrier of the 6.0 /im 
feature? Answering these questions will help us under- 
stand the nature of the HAeBe systems from which the 
PAH emission arises. 

Quantitative measures of the size and ionization frac- 
tion of the PAHs would be valuable in determining the 
degree of grain processing in the disks and the fraction 
of short-wavelength stellar photons that reach the outer 
disks. With the small size of the current sample, it is 
not possible to determine the origin of the correlations, 
although we expect that the full IRS sample of HAeBe 
stars will bring us closer to this goal. 



J. D. Bregman, D. M. Hudgins, and L. J. AUaman- 
dola provided useful comments as we prepared this 
manuscript. This work is based on observations made 
with the Spitzer Space Telescope, which is operated by 
the Jet Propulsion Laboratory, California Institute of 
Technology under NASA contract 1407. Support for this 
work was provided by NASA through contract number 
1257184 issued by JPL/Caltech. This research has made 
use of the SIMBAD database operated at the Centre de 
Donnees astronomiques de Strasbourg. 
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